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t> 

CO 
CN . 

The results and main steps of an analytic calculation of radiative-recoil correc- 
tions of order a(Za) 5 (m/M)m to the Lamb shift in hydrogen are presented. 
' The calculations are performed in the infrared safe Yennie gauge. The discrep- 

ancy between two previous numerical calculations of these corrections existing 
in the literature is resolved. Our new result eliminates the largest source of the 
D , theoretical uncertainty in the magnitude of the deuterium-hydrogen isotope 



shift. 



^ ■ I. INTRODUCTION 



The spectacular experimental progress in precise measurements of the energy levels in 
light hydrogenlike atoms achieved in recent years was matched by an equally impressive 
theoretical developments (see, e.g., the review in and references therein). Still, there is a 
number of unsolved theoretical problems, not least among them the magnitude of radiative- 
recoil corrections to the Lamb shift of order a(Za) 5 (m/M)m. This is the first nontrivial 
radiative-recoil correction, and it is generated by the radiative insertions in the exchanged 
photon lines and in the electron line. The correction generated by the one-loop polarization 
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insertions in the exchanged photon lines was independently calculated analytically in || 
and ||. The results of these works are in excellent agreement. Radiative-recoil corrections 
generated by radiative insertions in the electron line were obtained numerically in 
and in ||. The results of these works contradict each other. In this work we present an 
analytic calculation of the radiative-recoil corrections of order a(Za) 5 (m/M)m, and resolve 
the discrepancy between different theoretical results. 

Calculation of the contributions of order a(Za) 5 m to the Lamb shift is greatly facilitated 
by the applicability of the scattering approximation. These corrections are generated by the 
diagrams with at least two photon exchanges in Fig. [T], and naively one could expect that 
the diagrams with larger number of exchanges are also relevant. However, this does not 
happen. First, one has to realize that for high exchanged momenta expansion in Za is valid, 
and addition of any extra exchanged photon always produces an extra power of Za. Hence, 
in the high-momentum region only diagrams with two exchanged photons are relevant. 
Contribution of the low-momentum region is suppressed because the infrared behavior of 
any radiatively corrected Feynman diagram (or more accurately any gauge invariant sum of 
Feynman diagrams) is milder than the behavior of the skeleton diagram. Hence, unlike the 
leading contribution to the Lamb shift, diagrams with higher number of photon exchanges 
do not contribute to the corrections of order a(Za) 5 m, and it is sufficient to calculate only 
the contributions of the diagrams in Fig. |l] in the scattering approximation (for more detail, 
see, e.g., 0). 



diagrams with 
crossed exchanged 
photon lines 



FIG. 1. Electron-line radiative-recoil corrections 



In the scattering approximation contribution to the energy shift generated by the dia- 
grams in Fig. [I] is given by the integral 



AE = - t 1, 

7rn d 



/^4i Tf -[( 1+ 7o)V] -Tr[{l + 1 ,)H, V ] S l0 , (1) 



where m and m r are the electron and reduced masses, respectively, k is the momentum of the 
exchanged photon, and H^ v are the electron and the proton factors, respectively, and the 
Kronecker symbol Siq reminds that the radiative-recoil corrections of order a(Za) 5 (m/M)m 
are different from zero only for the S-states. 

The electron factor is equal to the sum of the self-energy, vertex, and spanning photon 
insertions in the electron line 

= + 2L^ U + L~ v , (2) 
and the heavy-line factor is given by the expression 

P + k + M P-k + M 
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where P = (M, 0) is the momentum of the proton. 

The expression for the energy shift in Eq.fll]) contains both recoil and nonrecoil contri- 
butions of order a(Za) 5 m. Nonrecoil correction of this order is well known from the early 
days of quantum electrodynamics, and in order to simplify further calculations we would 
like to obtain an expression for the energy shift free of such contributions. Let us notice to 
this end that the characteristic integration momenta in Eq.fll]) are of order of the electron 
mass, the lower momenta are suppressed by the radiative insertions in electron line, and the 
momenta of order of the heavy mass are suppressed by the high power of the integration 
momenta in the denominators. Suppression of the high integration momenta means that 
the radiative-recoil correction of order a(Za) 5 (m/M)m does not contain logarithms of the 
mass ratio ln(m/M) which could originate only from the wide integration region between the 
electron and proton masses m <C k <C M. Then we can kill all nonrecoil contributions and 
extract the first order in the mass ratio contribution simply by differentiating the integral 
in Eq.fll]) over the heavy mass M, and letting this mass to go to infinity afterwards. The 
integral in Eq.fll]) contains the heavy mass only in the heavy particle factor, and to extract 
the linear in the mass ratio term we make in the integrand the substitution 

irr[(l+ 7 „)ff„„] — 



P-k + M _ P + k + M 

7/i k 2 - 2Mk + i0 lu + lv k 2 + 2Mk + iO 1 * 1 



I A; 2 + — 
— ► - tt { k 2 g^g uQ - k ( g^K + g^k^) + k 2 g^ } — -j . (4) 

Due to the explicit factor 1/M before the braces, it is sufficient to consider the last factor 
in this expression only in limit k/M — > (we remind that the characteristic integration 
momenta are much less than the proton mass k <C M). Then the coefficients before the 
tensor structures in the last line in Eq.(|j) simplify dramatically. 

We are going to calculate the Feynman integrals in the polar coordinates in the four- 
dimensional euclidean space, where the vector k is parametrized in the form k Q = k cos 9 and 
|k| = ksmO. It is easy to see that the integrals with the coefficient before the first tensor 
structure in the braces in Eq.(|) should be calculated as principal value integrals over the 
polar angle 

U 2 y k l + im y cos 2 e-e 2 ( i x 

i^o (^0-4^2) (cos 2 ^ + e l Y ^co^Q 1 

where e = k/2M and p is the principal value symbol. 

The second term in the braces in Eq.(|]) is odd in ko, and it will give nonzero contribution 
to the integral over k only multiplied by another odd in ko term from the electron factor. 
But this means that the electron factor would effectively supply an extra power of ko in the 
numerator, and we can safely take the limit e — > in the coefficient before the second tensor 
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structure in Eq.(Q). The factor before the g muv term in Eq.(f4|) admits a smooth limit for 
e — > 0. Then the heavy particle factor turns into 



-Tr 



1 + 7o)#^] -> - — { fc 2 S^oS^o p(p 



(6) 



and the general expression for all radiative-recoil corrections of order a(Za) 5 (m/M)m ac- 
quires the form 



AE 



(Za) 5 ml r d 4 k 1 



7m 3 M J m 2 k A 4 



-Tr{ (l + 7 o) + 2L* + 



{ k 2 g^o9uo p{to) ~ ( 9^oK + 9uok^) jr + 9p» } 5 io ■ (7) 

This expression is much more convenient for calculations than Eq.(|l]) because it depends on 
the heavy mass only through the explicit factor m/M before the integral. 

The expression for the energy shift in Eq. (0) is linearly infrared divergent like I/7 where 
7 is an auxiliary infrared cutoff in integration over k. This linear infrared divergence is the 
price we have to pay for the simplicity of the scattering approximation. The point is that 
the expressions for the energy shifts in Eq.(|T]) and Eq.(^) contain not only corrections of 
relative order {Za) 5 but also the corrections of the previous order in Za. If we would not 
ignore small virtualities of the external electron lines and the external wave functions the 
naive infrared divergence would be regularized at the characteristic atomic scale 7 ~ mZa. 
We are not interested in the contributions of the previous order in Za, and will simply 
throw away linearly infrared divergent contributions in our calculations. Remaining infrared 
finite contribution is just the radiative-recoil correction of relative order a(Za) 5 . This 
strategy works because the radiative-recoil corrections under consideration do not contain 
logarithms of Za. However, individual diagrams could contain logarithms of the infrared 
cutoff which should cancel in the final result. This cancellation serves as an additional test 
of the correctness of all calculations. In the scattering approximation the diagrams in Fig. [I] 
form a complete gauge invariant set, and we can use arbitrary gauge for their calculation. 
In order to improve the low-momentum behavior of individual diagrams we use the Yennie 
gauge for the radiative photons. 



II. MASS OPERATOR CONTRIBUTION 

Let us consider first the contribution to the radiative-recoil correction of order 
a(Za) 5 (m/M)m generated by the diagrams with the self-energy insertions in the electron 
line in Fig. [I]. The renormalized mass operator in the Yennie gauge has the form (see, e.g., 



s » = r G 5 -™) 2 ft* 2 M( 2 3px 2 , n — \ • (8) 

Aix Jo m z x + \m A — p z ){l — x) 
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According to Eq.(|7]) respective contribution to the Lamb shift may be written as 



AEs 



3 a(Za) 5 mf 
A n 2 rfi M Jo 



dx 



d A k 



x 



in 2 k 4 A 1 A 



Tr 



(1 +7o) 7m ( P~ k ) lv 



{ & 2 g^agvo - ( fl^f + ft*)^) + £V } 



3 m, 

4 7T% 3 M 



5 ^ /'dx 
o 



x 



in 2 k 4 A 1 



Am + 2k Q + A; 2 — + m/c 2 1 , (9) 

k \k.zs J 



where 



Ai = m 2 x + 2p&(l - x) - k 2 {\ 



x 



x) 



-k 2 + 2mk + a\) 



and a 2 = m 2 x/(l — x). In the formulae below we will often use dimensionless mo- 
menta measured in terms of the electron mass k — > mk, and a shorthand notation 
C = a(Za) 5 /(ir 2 n 3 )(m/M)(m r /m) 3 m for the common normalization factor. 

As an example of our calculations let us consider evaluation of the contribution to the 
energy shift generated by the last most infrared singular term in the square brackets in 
Eq.® 



AE, 



3C 

~a 



dx 



d 4 k 



x 



m 2 k 2 (-k 2 + 2k + a 2 ) P ^k 2 



3C r 1 2 dk 2 2 sin 2 9 (k 2 + a 2 ) 

— ^ dx a x — - ^ dO ^ 2 + a 2y + 4k 2 cos 2 q vy^pg. 



(10) 



The integration in the last line goes over the four- dimensional euclidean space, and we have 
introduced an auxiliary infrared cutoff 7. Using the identity 



k 2 + a\ 



Ak 2 cos 2 6 



(k 2 + aj) 2 + Ak 2 cos 2 6 k 2 + ai (k 2 + a 2 ) [ (k 2 + a 2 ) 2 + Ak 2 cos 2 6 
we represent the integral in Eq . (|i~0|) in the form 



AE, 



3C /-i 2 f°° dk 2 2 r* Ja >5 
— / dx a x I — — au sm 
4 Jo J-y k % Jo 



k 2 + a\ P ^cos 2 6 



(11) 



4A; 2 



(k 2 + a 2 )[ (k 2 + a 2 ) 2 + Ak 2 cos 2 9 



where the principal value integral contains only a trivial dependence on the angles. 
All principal value integrals we need in this work, in particular the integrals 



2 r n / 1 

- / dO sin 2 6p(—- 

7T JO x -COS z O 



- 2 



(12) 



9 rn , 1 

- / d9 sm 4 6p( 

7T JO V COS 

may be obtained from the basic integral 

r d6 p( 



2 



3 , 



cos 2 6 







with the help of algebraic transformations. 

Now we can easily complete calculation of the integral in Eq.([TD|) 



3C yi 



rfx <! 2 In 4 - 2 In 

7 



1 + of 4 1 

5 — + — arctan — 

af a d\ 



C{ 3 In 



7 



(13) 



3tt 2 



Other integrals in Eq.(^) are calculated in the same way, and we obtain the total self- 
energy contribution to the energy shift in the form 



5E, 



7T 



/m r \ 3 




— to 


9 ln- 


\m J 


L 7 J 



(14) 



III. VERTEX CONTRIBUTION 



We calculate the contribution to the Lamb shift generated by the vertex insertion in the 
electron line in Fig. [I] with the help of the compact expression for the electron-photon vertex 
in the Yennie gauge used in our earlier work on the radiative corrections of order a 2 (Za) 5 m 



n=! A' 



(15) 



where 



F« = 37 M f k 2 - 2p k + (2 - x) A 



— x 



37„ (p + to) 7 M (p - k + to) 7" 



6(p - k) Q 7 M + 7 a Q 7 M (p -k + to) 7" + 7 Q (p + m) 7^ Q 7° 



+ 2 7^ (p - k + to) Q + 2 Q (p + m) 7^ ] + ( 2 Q 7/t Q + Q 2 7m 



(16) 



F< 2) = 2x(l - x) 7^ (p - fc + to) f Q p Q - p Q 



(17) 



A = m 2 x + 2p/c(l — x)z — A; 2 2;(l — rrz) = z(l — xz)( a 2 + 2kob — A; 2 ) 



(18) 



and Q = —p + kz, p 2 = m 2 . 

According to Eq.(^) radiative-recoil contribution of order a(Za) 5 (m/M)m generated by 
the vertex insertion has the form (we use dimensionless integration momenta below) 
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AE, 



a(Za) 5 m /m r \ 3 r 1 , r 1 , /■ d 4 £; 1 I 
2 „ „ — — I — I m ax dz 



7r 2 n 3 M \ m 



o Jo 



tir 2 k* k 2 - 2k 



n=l 



where 



(19) 



K (n) = -Tr 



(l+7o) (p-k + 1) 



Vl n) = - a Tt (1 + 7 o)fW (p-k + 1) j v ( - <^ - g vQ k 



M 7 ' 



K, (n) = -Tr 
2 4 



(1 + 7o ) F< n > (p - k + 1) lv ] k 2 g, g u0 . 



Calculating traces and contracting the Lorentz indices we obtain the numerator factors in 
the square brackets in Eq.(|T9|) 



V n (1) + V> 



Ax(l - x)k 2 



( 1 

fc 



+ 2x(5 - 9z + 6z z + Axz - Axz 2 )k 2 k + f- 



fc J 



+ (A; 2 -2A;o){ [ - 6(1 - 2z) + 2z(-l - 10 z + 2x + 2xz) 



+ 



6(1 - 2z) + 2x(l + z + 6z 2 + 2xz - Axz 2 ) 1 fc 2 pf i) 



(20) 



+ 



6(1 - 2z) + 2x(-5 + 10^ - Axz) k + 3(1 - 2z) + x(z + Qz 2 - Axz 2 ) — \ 



4x(l - x)2 k 2 { 4(1 - z) 



1 + k 2(1 - z) (2fc 2 - ^) - 



+ (A; 2 - 2k ) [ _ 22 + 2(1 - 22) — + zArW-^) 1 } 



(21) 
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As a simple example let us consider calculation of the contribution to the energy shift 
generated by the terms in the last line in Eq. (f20D 



1 r 1 r d 4 k 1 



6(1 - 2z) + 2x(-5 + lOz - Axz) k 



+ 



3(1 - 2z) + x(z + 6z 2 - 4xz 2 



k 



2C C dx f 1 dz U 



dk 2 - / de 



o Jo z(l — xz) Jo n Jo (k 2 + a 2 ) 2 + Ab 2 k 2 cos 2 9 



sin 2 9 



| [6(1 -2z)+2x(-5 + IQz-Axz)] cos 2 9 + [3(1 - 2z) + x(z + Qz 2 - Axz 2 )] | . (22) 

The last integral as all other integrals in this paper may be written as a linear combination 
of the integrals of the form 



oo 9 ft 

dk 2 - / d9 

7T Jo [{k 2 + a 2 ) 2 + Ab 2 k 2 cos 2 9}p 



2\n 



sin 2 9 cos 2 ' 9 (k 2 ) m (k 2 + a 2 ) 



(23) 



where l,m,n = —1, 0, 1 and p = 1,2, 3. For I = —1 the integrals with cos 2 9 in the denomina- 
tor should be interpreted as principal value integrals. All these integrals may be calculated 
in terms of four standard functions of the parameters a, b 



a 2 + b 2 a 
L = In — , Li = 1 - — L 



o- 



a' b 2 
The result of integration for our example in Eq.(p2|) is 

1 



2a 2 T b b 

1 rr-^i> — arctan— . 

b l a a 



AE\ = 2C C dx t dz ■ 

A Jo Jo z(l - x) 



6(1 - 2z) + 2x(-5 + lOz - Axz)\ [-L 



(24) 



l - Ll 
2 



+ 



3(1 — 2z) + x(z + 6z 2 — Axz 2 ) (—L + — arctan - 



a(Za) 5 m ( m r \ 3 
2 \ . — m 



7r 2 n 3 M 



rn 



11 



C(3) + 7tt 2 In 2 



53tt 2 
12 



+ 7 



(25) 



In the contribution to the energy shift in Eq. (|22| ) the electron denominator canceled with 
a similar term in the numerator, and as a result all denominators were combined with the 
help of only two Feynman parameters x, z. Let us turn now to a simple example when such 
cancelation does not take place, and we need to introduce a third Feynman parameter t. 
The contribution generated by the first term in Eq . (|20|) may be written in the form 



S 



1 



^ = 8C lo^Io^I^ k 2 (k 2 -2k ) 



x(l — x) 
A 



-8C ! dx [ dz(l-x)a 2 [ dt 
Jo Jo Jo 



d_ 

da? Jo 



oo O 

dk 2 - 



k 2 + a 2 

d0sm 2 6-—: ^ 5-r, (26) 

77 ./ii (k 2 + a 2 ) 2 + Ab 2 k 2 cos 2 6 v 7 



where a 2 = a 2 t and fo t = 1 — z(l — z)a 2 t. 

The integral over angles and momenta in Eq.(^) is just of the standard form Eq.(p3|). 
and may be easily calculated. Integration over t is facilitated by the simple observation 
that the t- integral may be written as an integral over a 2 with the upper limit a 2 . After this 
transformation all dependence of the new integral on x is hidden in this upper integration 
limit, and we can get rid of the third Feynman parameter integrating by parts over x (for 
more details see [[?]) 



AE , 



8C f 1 dx f 1 dz(l -x)a 2 C 5 L 
Jo Jo 



b 2 



0/ 



8C 



\ dx \ dz(l-x) / L ot = AC / dx / — L = 2C 

Jo Jo Jo of Jo Jo z 



1L-2 
2 



(27) 



Contributions of the other terms in Eq.(EO) are calculated in the same fashion, and the 
total contribution of all terms with n — 1 in Eq.(|19|) is equal to 



AE 



(i) 



^ a(Za) 5 m /m r \ 3 



7r 2 n 3 M Km 



91n- + |C(3) 
7 8 



9tt 2 



In 2 



3tt 2 



+ 4 



(28) 



For the total contribution of all terms with n = 2 in Eq.(|19D we obtain 



AE 



.2, = ._ /'(^>)'"' "> ('">■)''' u 

7i 2 n 3 M \ m 



21 



C(3) 



13tt 2 



ln2 



97T 2 

~8 



35 
T 



(29) 



Total vertex insertion contribution to the radiative-recoil correction of order 
a(Za) 5 (m/M)m is given by the sum of the contributions in Eq.(P%D and Eq.(E^) 



AE, 



a(Za) 5 m f m 



71 



2 n 3 M \ m 



m 



18 In- + 6C(3) - 27T 2 ln2 



3tt 2 



7 



19 

T 



• (30) 



IV. SPANNING PHOTON CONTRIBUTION 

The contribution to the Lamb shift generated by the spanning photon insertion in the 
electron line in Fig. [I] is calculated with the help of an explicit expression for the jellyfish 
diagram. The small k behavior of the jellyfish diagram is one of our primary concerns in 
further calculations, since the contribution of the previous order is connected just with this 
infrared region. The jellyfish diagram is finite at small k in the Yennie gauge, and this is 
one of the reasons why we are working in this gauge. We need such representation for the 
jellyfish diagrams where not only the diagram as a whole, but all entries are finite at k — 0. 
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The compact expression for the jellyfish diagram with such properties was used in our earlier 
work on the radiative corrections of order a 2 (Za) 5 m || 



r s = — 



I dx I dzxil — z) 
Jo Jo n=1 



3 M (n) 

~^ fJbV 



(31) 



where 



M (1) = -2N {a \ 



(32) 



M$ = NjS + 2(1 - x)jv£> + 



/IV 1 [W 



Mj$ = -2(1 -x)Nj$, 



/IV 



NjS = 2m[Q 7/ ,(p + xQ - k + m) 7l/ + 7m (j3 + a;Q - k + m)j u Q] 



(33) 



-12(1 - 2x)m 2 ^^Q^ u - 2xQj u (p + xQ - k)^^Q + AxmQ 2 g llu 



+ [2xQ 2 + 8(pQ)h»{p + xQ-k + m)^ u , 



N$ = 4(pQ) 7m P7„ + 2m 2 7 M Q 7 „ 
iVW = 8[( P Q) 2 - m 2 Q 2 }^(p + xQ-k + m) 7 „ 



and a 2 , 6 and Q where defined in Eq. (|i~8|) . Notice that all Mf$ are infrared finite even at 
k = 0. 

General expression for the energy shift induced by the spanning photon insertion has the 
form (see Eq.(f7])) 



a(Za) 5 m ( m r 
7r 2 n 3 M \ m 



m dx dzx(l — z) — - 
Jo jo J in z 



d 4 k 1 

2 ~ui 



n=\ 



T, 



(n) 



+ T i n) F + T 2 n) p(^) 
k v fco y 



(34) 



where 
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T 



(n) 



- Tr 
4 



(1 + 7o) 



9nv 



n 



- Tr 
4 



(1 + 7 o) MfJ 1 ( - g^K - 9v0 k^ ) 



(i i i 



T 



(n) 



- Tr 
4 



(l + 7b)MW A; 2 



(35) 



Calculating traces and contracting the Lorentz indices we obtain the numerator factors in 
the square brackets in Eq. (|34]) 



71(1) I t (± > — + T w <ni — 

h + 11 h + 12 P U 2 



,(1) 



,(1) 



24 (1 - x) +4 (-3 + 2xz) k + 2 (-3 + 2xz) — 

k 



+ 6 (-3 + 2x) k 2 pija) 

^0 



(36) 



T (2) + T| 



i2) — + T (2) nf— 
k + 2 P U 2 



12x - 8z(l - x)(4 - x) k + %z(2 - x - xz) kl 



+ Az(A + x - 4x2 + 2x 2 z) k 2 + 8z(l - x) (4 - x) — 

fe 



+ 6xA; 2 p(-V) + 2z(2 + 2x + 3xz - Ax 2 z) k 4 p(-^ 



k\> 



(37) 



t (3) + I + rf 



- 16(l-x)z 2 k 2 2(-l + 2x) + 2(l-zz)Jfe 



+ (1-xz)^- + (2-x) fc 2 p(-^) 
fc v 



(38) 



We will use the integral 
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. / „ r 1 , r 1 , . , f d A k k 2 — 2xz / 1 \ , , 

AE 3 = l 6C ^/ o ^(l-z)/-- — (39) 

describing one of the contributions with n = 3 in order to illustrate one more subtlety 
encountered in our calculations. This infrared divergent integral, contains in the integrand 
the term x/A 3 . At small k — > the denominator A — > x, and integration over x becomes 
too singular. The singular factor pil/k^) in the integrand makes things even worse, and 
we risk to end up with a divergent integral over the Feynman parameter x, instead of an 
infrared divergent integral over the momentum. While linearly infrared divergent integrals 
over k have a transparent physical interpretation as contributions of the previous order in 
Za, divergent integrals over x could lead to uncontrollable contributions. We separate the 
infrared divergent part of the momentum integral with the help of the identity 

yi x 1 2z z 2 (l-z)k 2 

Z Jo X A 3 ~ k 2 - 2k + 1 - k 2 z(l - z) + [1 - k 2 z{l - z)] 2 1 () ) 



2x\ 

a^J 



Then the singular integration over momentum decouples, and we easily obtain 

A£~ = 16C 



4 3 
37 ~ 4 



(41) 



After tedious calculations we obtain 



AE (i) = a{Zaf m /m r y 
7r 2 n 3 M V m 



m 



9 In 



7 



^C(3) + -77- In 2 



15vr 2 159 
+ T 



(42) 



a(Za) 5 m ( m r \ 3 



7r 2 n 3 M 



rn 



m 



12C(3) 



3tt 2 



- 16 



(43) 



AE 



(3) 



a(Za) 5 m /m 



vr 2 n 3 M 



r?) 



75 
T 



3vr 2 n 9tt 2 

C 3 In 2 + 

SV ; 2 16 



113 



(44) 



for the contributions with n = 1,2,3 in Eq. (p4[) , respectively. Let us emphasize once again 
that we have thrown away linearly infrared divergent term I/7 in AE E , which corresponds 
to the contribution of the previos order, but preserved all logarithmically divergent terms 
which should cancel automatically in the final result for the energy shift. 

Total contribution to the energy shift generated by the spanning photon insertion in 
Fig. [I] is equal to 



5E^ 



aiZaf m (m 



n 2 n 3 



M 



rn 



m 



, 1 9 
9 In 

7 2 



(45) 
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V. SUMMARY 



Collecting all contributions to the energy shift in Eq.(|TJ|), Eq.(^), and Eq. we see 
that all logarithmically infrared divergent contributions cancel in the sum, and obtain the 
total radiative-recoil correction of order a(Za) 5 (m/M)m 

. / . o, 37r 2 \ a(Za) 5 m /m r \ 3 , , 

AE = 6C(3)-2 I »ln2 + — -14 -L^ - (_) m (46) 



a(Za) 5 m /m r \ 3 
« -13.0676322 . . . \ / — I — m. 

7i n M \m J 

This result is in excellent agreement with the numerical result in |3[], and this resolves 
the long standing discrepancy on the magnitude of the radiative-recoil corrections of order 
a(Za) 5 (m/M)m to the Lamb shift. When this paper was in preparation we learned that 
the same result was just obtained in the NRQED framework by Czarnecki and Melnikov ||. 
Numerically, the correction in Eq.(f46j) contributes 



AE(1S) = - 13.43 kHz (47) 

to the IS* Lamb shift in the ground state of hydrogen. The discrepancy between the the- 
oretical predictions for the IS Lamb shift calculated according to and || is about 6 
kHz. This discrepancy is not too important for the IS Lamb shift measurements, since the 
error bars of even the best current experimental results are still a few times larger (see, e.g., 
review in 0). What is much more important from the phenomenological point of view, the 
radiative-recoil correction is linear in the electron-nucleus mass ratio, and it directly con- 
tributes to the hydrogen-deuterium isotope shift. The discrepancy between the theoretical 
values of the isotope shift calculated according to |§-§ and || is about 18 times larger than 



the experimental uncertainty 0.15 kHz of the isotope shift | ID| . Thus the analytic result in 
Eq.(f46"D eliminates the largest source of the theoretical uncertainty in the magnitude of the 
deuterium-hydrogen isotope shift. One can use the new radiative-recoil correction, and the 
latest experimental data in order to obtain a new value for the difference of charge radii 
squared of the deuteron and proton, but we will not enter in the detailed discussion of the 
phenomenological implications here, since they were exhaustively discussed in our recent 
review lllfl. 



ACKNOWLEDGMENTS 

We are grateful to A. Czarnecki and K. Melnikov for communicating us their results 
prior to publication. 

This work was supported by the NSF grant PHY-0049059. 



13 



REFERENCES 



[1] M. I. Eides, H. Grotch, and V. A. Shelyuto, preprint PSU/TH/226, |hep-ph/0002T58| , 

February 2000, Physics Reports, in print. 
[2] M. I. Eides and H. Grotch, Phys. Rev. A52 (1995) 1757. 
[3] K. Pachucki, Phys. Rev. A 52, 1079 (1995). 
[4] G. Bhatt and H. Grotch, Phys. Rev. A 31, 2794 (1985). 
[5] G. Bhatt and H. Grotch, Phys. Rev. Lett. 58, 471(1987). 
[6] G. Bhatt and H. Grotch, Ann. Phys. (NY) 178, 1 (1987). 

[7] M. I. Eides, S. G. Karshenboim, and V. A. Shelyuto, Ann. Phys. (NY) 205 (1991) 231. 
[8] M. I. Eides and V. A. Shelyuto, Pis'ma Zh. Eksp. Teor. Fiz. 61 (1995) 465 [JETP 

Letters 61 (1995) 478]; Phys. Rev. A52 (1995) 954. 
[9] A. Czarnecki and K. Melnikov, preprint SLAC-PUB-8730, |hep-ph/0012053; , December 

2000. 

[10] A. Huber, Th. Udem, B. Gross et al, Phys. Rev. Lett. 80 468 (1998). 



14 



